Bond-Ordering Model for the Glass Transition 
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We consider the idea of bond ordering as a model for the glass transition; a generic covalently bonded 
liquid may substantially reduce its energy through bond ordering, without undergoing significant 
structural order. This concept is developed for a model system using quantities such as a bond order- 
parameter and susceptibility which provide a new identification for the calorimetric glass transition 
temperature. Monte Carlo simulation results exhibit bond ordering at intermediate temperatures 
uncorrelated with any long-range structural ordering. Also discussed are various other implications 
of bond-ordering model for the glass transition. 



I. INTRODUCTION 

The problem of glass formation is a classic one in 
physics and materials science. Viewpoints break down 
into kinetic explanations of an arrested liquid, and phase 
transitions with some kind of order parameter. Discus- 
sions have approached the glass transition in terms of 
the fragmentation of small structural units,tl the agglom- 
eration of clusters through chemical bonds ,□ or corre- 
lations between different metastable equilibrium states 
representing distinct configurations or rearrangements 
of the systemE In addition, attempts have been made 
to incorporate frustration in local bond ordering in a 
glass through the introduction of a local order param- 
eter describing locally preferred arrangements of liquid 
molecules, which, in general, are not consistent with t 
crystallographic symmetry favored by density ordering 
In this model the frustration arises from competition be- 
tween density ordering and local bond ordering, explain- 
ing why some molecules crystallize easily without vitri- 
fication, while others easily form glasses without crys- 
tallization. Another approach, using computer simula- 
tions, introduces a displacement-displacement correla- 
tion function as a measure of the local configurational 
order which grows as the temperature is lowered toward 
the glass transition.^! Finally, the role of fluctuations has 
been reviewed in general fashion with both theoretical 
and experimental evidence for heterogeneity at the glass 
transition^ We would like to introduce a model of po- 
sitional glass which incorporates many of these ideas in 
simple form. 

First we recall some of the basic structural features of 
typical network-glass-formers like silica. Vitreous silica 
is made up of SiC>4 tetrahedra just like the crystalline 
phase, quartz. Thus, molecular structure for the two dif- 
ferent phases of silica, i.e. crystalline and amorphous, are 
very much alike, with the major exception of n.n.n. dis- 
tance Sii-Si2 which may be specified in terms of the Si-O- 
Si bond angle. The distribution of Si-O-Si bond angles 
(3 as determined by Mozzi and Warren for vitreous sil- 
ica is shown in Fig. 1.0 From the bond angle distribution 



V(/3), it is clear that the Sii-Si2 distances vary signifi- 
cantly, which in fact is the primary source of topological 
disorder. Yet the bond angle (3 is confined to within 
roughly 30° of its most probable value viz 144°. Clearly, 
the variation in (3 is large enough to suppress the long- 
range order (LRO) characterizing the crystalline phase, 
and yet small enough to maintain the medium-range or- 
der (MRO) which extends to about 10RA. Hence, we 
note that ordering in n.n. and n.n.n. distances or bond 
lengths is well maintained in spite of the fact that vit- 
reous silica is devoid of any substantial structural or- 
der. Similar structural properties are observed in various 
other network glasses such as B2O3 and GeC>2. This ob- 
servation is striking enough to indicate a strong role for 
the ordering of bonds in a proper microscopic model for 
the glass transition. The bond ordering can be viewed as 
being brought about by the local reorientations of molec- 
ular clusters as suggested in the theory of Adam and 
GibbsE The idea is that supercooled liquids do not nec- 
essarily need to undergo structural ordering in order to 
achieve local equilibrium. In fact, local rearrangements 
of molecular clusters can lead to substantial lowering in 
the internal energy of an entire sample through reduc- 
ing bond energies at the local level, which is the primary 
reason for the ordering of n.n. and n.n.n. distances, or 
alternatively the bond-angle degrees of freedom, in amor- 
phous materials. 

A model for the glass transition incorporating small 
structural units or fragments was proposed by M. Suzuki 
et aln. In this fragmentation model noncrystalline solids 
are assumed to be assemblies of pseudo-molecules — a 
pseudo-molecule being a cluster of atoms having a disor- 
dered lattice in which there are no definite defects such 
as under- or over-coordinated atoms. As temperature 
increases, bond breaking intensifies at the boundaries of 
such clusters where the bonds tend to be weak. The bond 
breaking mechanism arising from the thermal excitation 
of electrons from bonding to antibonding energy states, 
causes the noncrystalline solid to fragmentize, with the 
average fragment size decreasing as the temperature in- 
creases. Consequently, material begins to show viscous 
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flow when average fragment size reaches a critical value. 
The model is shown to have some success in describing 
the temperature dependence of viscosity and the varia- 
tion in transition temperature with heating rate for a-Si. 
The origin of pseudo-molecules in the cooling process, 
however, is not addressed in the fragmentation model. 

In the following we consider a microscopic model for 
the glass transition which is general enough to be ap- 
plicable to various types of glass-forming systems with 
widely differing bonding schemes and chemical composi- 
tions. We would like to treat glass transition as a phe- 
nomenon that finds similar description whether a liquid 
is cooled or heated through the transition. The focus 
of bond-ordering model is the bonds linking neighboring 
atoms rather than the atoms themselves. In other words, 
a bond is treated as a distinct object in its own right, pos- 
sessing internal degrees of freedom or electronic states. 
The internal state of a bond is governed by the separa- 
tion of the participating atoms. The term bond-ordering 
refers to the process of relaxation of bonds into their low- 
lying internal energy states, facilitated by the cooperative 
rearrangements within molecular groups. Bond ordering, 
therefore, may be viewed as some form of ordering in en- 
ergy space. The important point that we would like to 
bring home, however, is that such an ordering can be 
achieved without need for any significant structural or- 
der. To this effect, we have provided results from Monte 
Carlo (MC) simulations of the model Hamiltonian which 
couples the coordinates of ions to the electronic states 
of electrons, as for a typical covalently bonded network 
glass. The simulations make clear the possibility of local 
ordering of bonds, uncorrelated with any kind of long- 
range structural ordering. 

Sec. II of this article contains the theoretical back- 
ground regarding the covalently-bonded systems. In 
Sec. Ill we describe the model Hamiltonian. Some new 
definitions in terms of mathematical expressions for vari- 
ous structural and bond-related quantities of interest, are 
introduce in Sec. IV. The simulation procedure and re- 
sults are discussed extensively in Sec. V. Sec. VI contains 
the concluding remarks and a summary of the main ideas 
introduced in this article. 



II. THEORETICAL BACKGROUND 



which describe the interactions between atoms. The oc- 
cupation of the bonding state by one to two electrons 
results in a lowering of the electronic energy relative to 
the atomic levels. This covalent bond provides the struc- 
tural constraint on the relative position of the two atoms. 

In order to be more specific we consider a general 
Hamiltonian for a covalently bonded network glass. We 
start first with a general formulation of the Hamiltonian 
for a typical glass-forming liquid. Although these sys- 
tems are often treated in a classical formulation, in fact, 
the bonding between atoms in a network glass, as well as 
the coordination in metallic glasses arise from quantum 
mechanical considerations. The many-body Hamiltonian 
may be written in the form: 

H = Ha + Hi e + H ee , (1) 

where Ha is the interaction between ions in the liquid, 
Hi e is the ion-electron interaction, and H ee is the inter- 
action between electrons. The usual approach to solving 
this type of many-body problem is to exploit the separa- 
tion of time scales inherent in the Born-Oppcnheimer ap- 
proximation: solve for the electronic states regarding the 
ionic coordinates as parameters, then vary the ionic co- 
ordinates to minimize the energy or extract the electron- 
phonon interaction, etc. For a covalent material like most 
network glasses, the electronic states are reasonably ap- 
proximated by forming bonding and antibonding states 
from the atomic s and p orbitals. The available elec- 
trons are then apportioned to the lowest bonding states 
to obtain the ground state of the system. The bonding 
and antibonding energies, depend on distance or separa- 
tion as shown schematically in Fig. 1: the bonding state 
has a minimum at an ion-ion distance r , while the anti- 
bonding state is repulsive at all distances. This may be 
thought of as a tight-binding approximation to the ac- 
tual electronic structure of a glass. A bond, is a bonding 
state occupied by two electrons and is strongest (or has 
lowest energy) when the distance between ions is near 
tq. On the other hand, a bond missing an electron due 
to a thermal fluctuation or a transition of the electron 
to the antibonding state, corresponds to a broken bond. 
In this tight-binding representation of bonds, the Hamil- 
tonian may be expressed as in Eq. (1), a more detailed 
expression of which is derived in the next section. 



Glass-forming liquids include covalently-bonded net- 
work glasses, coordination-constrained metallic glasses 
and systems made of more complex organic molecules 
or polymer. Here we take the view that the interaction 
between constituent atoms or molecular units is given by 
an effective potential which characterizes their positional 
and orientational interaction. In particular for covalent 
materials the interactions may be reasonably described 
in terms of basis states of degenerate s and p orbitals on 
each atom. Linear combinations of orbitals on neighbor- 
ing atoms then lead to bonding and antibonding states 



III. MODEL HAMILTONIAN 

The potential energy of a pair interaction may be ex- 
panded in the displacement x from the equilibrium: 



U{x) = U Q + -US 



3! 



(2) 



where the coefficients are evaluated at equilibrium sepa- 
ration. Dropping the constant term which plays no role, 
and in the harmonic approximation, we get: 
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U(x) 



-u" 



(3) 



For our lattice model we identify x = | i?j — Rj | , where Ri 
and Rj arc the displacements from the equilibrium of the 
atoms assigned to the zth and jth n.n. sites, respectively. 
Hence, with the idealizing constraint that all displace- 



ments are of same magnitude \R\ 
we have: 



\R, 



\R4 



U (\Ri - Rj\) = \R\ 2 U£ (l-Ri-R 



\R\, 



(4) 



Letting J = \R\ 2 Uq, we obtain the following Hamil- 
tonian for a system of atoms interacting through n.n. 
coupling: 



H = -J 



<i,j> 



(Ri ■ Rj) 



(5) 



This is clearly of the same mathematical form as the 
g-state clock model Hamiltonian, when the displacement 
degrees of freedom are taken to be discrete. 

The interaction Hamiltonian we used for our MC sim- 
ulations, is the following: 



n = - j • R j + 1 

<i,j> 

= —j (cos% + 1) Uij , 



(6) 



<h3> 



where 6ij = =g-(si — Sj), with s,, Sj = 1,2, ... , 8. The Sj's 
are integer labels for various possible displacements of 
an atom from its equilibrium lattice site i . The quantity 
riij should be regarded as the bonding-electron occupa- 
tion number for a (possible) bond, linking n.n. sites i 
and j. A bond may or may not be broken depending 
on whether the corresponding takes on the values 
0, or 1, respectively. The value taken by riy depends 
on the number of bonding-electrons made available to 
the system — this is an input parameter to the simulation 
code — and the relative values of energies of the n.n. in- 
teractions. Bonds with lower values of energy are more 
likely to have bonding-electrons. Interesting effects are 
observed with a bonding-electron (hole) concentration of 
about 60% (40%) and that is what we report later in this 
article. As we shall see shortly, the net effect of holes is 
to suppress the AF-like transition occurring for clock 
models with q > 4, and hence allowing us to observe 
the behavior of a disordered system with lowering of the 
temperature. 



A. Bond Order-Parameter 

We introduce bond order parameter (or bond magne- 
tization) as a measurable physical property that indi- 
cates the extent of bond ordering prevailing in a physical 
system. We begin with considering a two dimensional 
system of 4-fold coordinated atoms interacting through 
n.n. coupling of strength J. The bond magnetization of 
such a system is significant when bonds are in their low- 
lying energy states (as for a bond-ordered low tempera- 
ture phase), and negligible when bonds are distributed 
among all possible energy states with uniform probabil- 
ity which is indeed the case when the thermal energy is 
far in excess of the coupling strength J. To construct an 
expression for the bond magnetization of such a system, 
every n.n. pair of atoms is characterized by a vector By 
the purpose of which is to characterize the interaction 
energy. Bij is specified via the angle it makes with an 
arbitrary fixed axis in the plain of the system, (f>ij. The 
angle 4nj is written in terms of bond energy as: 



J 



-J < Eij < J , 



(7) 



where J is n.n. coupling strength, is the bond energy 
and — 7r < (fiij < 7T . An expression that fulfills all the 
requirements of an extensive bond magnetization is the 
following: 




2N 2N 
<i,j> <m,n> 



1/2 



(8) 



where angular braces stand for thermal average and N 
is the total number of atoms in the system. In terms of 
bonds' energies, the expression for intensive bond mag- 
netization of the system; m\, = Mb/2N , may be written 
as: 



1 

2N 



2N 2N 



1/2 



Y ( cos j( e y ™ £m ™) ) 



(9) 



<i,j><m,n> 

The normalization is chosen such that < < 1 . 



IV. SOME NEW DEFINITIONS 



B. Bond Susceptibility 



In Sec. I we defined bond ordering as that process in- 
volving relaxation of bonds into their low-lying internal 
energy states. There are a few more physical parameters 
relevant to our discussion which are defined in this sec- 
tion and are considered in context of the glass transition 
phenomenon. 



Bond susceptibility is the response function associated 
with bond magnetization Mb and its thermodynamic con- 
jugate field Hb, called bond ordering fieldJa The exact 
physical nature of is not yet known to us; however, 
bond susceptibility can still be calculated without ap- 
pealing to a knowledge of Hi,. 
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The change in Gibbs free energy in an infinitesimal 
bond ordering process may be written in terms of the 
newly introduced parameters Mb and Hb, as in the fol- 
lowing: 



dG = -SdT-MhdH; 



b ■ 



(10) 



Eq. (|l^) can serve as starting point for incorporating Mb 
into the thermodynamics of disordered systems. Bond 
susceptibility apart from a normalization will therefore 
be given by: 



Xb 



8Mb 

dm 



d 2 G 
dHl 



(11) 



Starting with Eq. ([11]), one may readily obtain an expres- 
sion for the bond susceptibility in terms of fluctuations 
in the bond magnetization: 



Xb = ((\M b 



- (\M h 



/2Nk B T , 



(12) 



where Mb = Y^<i j> Bij ■ Eq. (|l^ ) is properly normal- 
ized to the number of n.n. bonds 2N for a system of N 
atoms with periodic boundary conditions, and is used for 
calculating bond susceptibility via MC simulations. 

This rather abstract entity, the bond susceptibility, de- 
scribes the tendency of a system for bond ordering and 
provides the basis for a new identification of glass tran- 
sition temperature. 



C. New Identification for Tg 

Following the previous discussion one can trace the ori- 
gins of MRO characteristic of the vitreous state in the 
local ordering of bonds, which becomes most intense at 
some particular temperature T g . This brings us to an- 
other identification for the calorimetric glass transition 
temperature, which particularly applies to fragile and in- 
termediate class of the glass-forming liquids £j In glass 
transition region, bond susceptibility of a supercooled 
glass-forming liquid reaches a maximum. This also im- 
plies that specific heat must display a maximum in the 
glass transition region because of the large energy fluc- 
tuations associated with the intense ordering of bonds. 
In view of this, the transition peak in the experimentally 
measured specific heat of various fragile and intermedi- 
ate glass-formers can be regarded as an artifact of in- 
tense bond ordering or strong fluctuations in the number 
of bonds in each energy state occurring at glass tran- 
sition. Hence, we propose a new identification for T g 
as that temperature corresponding to maximum of bond 
susceptibility of a supercooled liquid nearing configura- 
tional arrest. 

Having given an explanation for the specific heat peak 
at the glass transition, it doesn't seem improper to con- 
sider the unexpected linear behavior of specific heat at 
very low temperatures.lill Work by Hunklinger et al. gives 



strong evidence that the low temperature anomalous 
properties of amorphous materials arise mainly from two- 
level systems and not from the multi-level vibrational de- 
grees of freedom associated with the atoms.Ej The energy 
gaps A of the two-level systems are supposed to vary 
with uniform probability in some range < A < Ao, 
where Ao « lX.tLiSliJ These low energy excitations may 
be attributed to the bond ordering process at low tem- 
peratures. In disordered systems, whether supercooled 
liquids or glass, bonds continuously relax into more sta- 
ble internal energy states with lowering of the tempera- 
ture. At very low temperatures certain number of bonds 
may be seen to act like two-level systems with varying en- 
ergy gaps. Bond ordering process at low temperatures is 
therefore a possible explanation for the low temperature 
anomalous properties exhibited by amorphous materials. 



D. Structural Magnetization and Susceptibility 

Structural magnetization is a measure of the conven- 
tional LRO a system may possess. For the system we are 
considering the displacement of an atom from its des- 
ignated equilibrium lattice site i, is characterized with 
Ri. Further idealizing the system by assuming that dis- 
placements are of same magnitude \R\, we can express 
extensive structural magnetization as in the following: 



M. = 




(13) 



where Ri — Ri/\R\ . The intensive (or normalized) 
structural magnetization is given by; m s = M s /N, where 
< m s < 1 . The analogy between structural magneti- 
zation and magnetization of a magnetic system, is rather 
obvious. In fact, we can use this analogy to express struc- 
tural susceptibility in terms of fluctuations in structural 
magnetization, as follows: 



X. 



= {(\M s \ 2 )-(\M s \) 2 )/Nk B T 



(14) 



where M s = (J^-Rf , J2i^f)- Structural susceptibil- 
ity is the response function describing the tendency for 
structural ordering. 



E. Order-Parameter for Glass 

We would like to address at this point a possible or- 
der parameter (or order parameter density) for super- 
cooled liquids and glass, that also serves as yet another 
distinction between fragile and strong classes of the glass- 
forming liquids. This will simply be the bond magneti- 
zation if one is considering strictly a liquid <-> glass tran- 
sition. 

Many of the theories describing glass transition phe- 
nomenon, assume that there is a single parameter which 
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characterizes, glass. This assumption is believed to be 
inaccurate. c3 Prigogine and Defay have shown that in 
general the ratio of the discontinuities in second-order 
thermodynamic quantities; isothermal compressibility 
heat capacity at constant pressure, and coefficient of 
thermal expansion: 



R = 



Ak t AC p 
TV{Aa) 2 



(15) 



is equal to unity if a single order parameter characterizes 
the underlying thermodynamic transition, but if more 
than one order parameter is involved, then R > 1 .L£l 
The latter seems to describe most glasses. In view of 
this, we consider a two parameter description of glass 
which involves two of the parameters described earlier, 
namely, the structural magnetization m s and the bond 
order parameter rrib. We require for amorphous solids 
that the structural LRO vanishes while the bond magne- 
tization remains significantly large. The requirement of 
the vanishing of structural LRO, is meant to character- 
ize the liquid-like attributes of the amorphous systems. 
Yet large values for bond order parameter is a solid-like 
attribute that should serve to distinguish the glass from 
the liquid phase. 

Clearly, in case of strong glass-formers characterized 
with strong covalent bonds, the values for bond order 
parameter above and below the transition must be quite 
comparable. However, in case of fragile systems the liq- 
uid undergoes substantial bond ordering at the transition 
mainly due to the nondirectional nature of their chemi- 
cal bonds. As a result bond magnetization is expected to 
vary rather significantly for the fragile class. Originally, 
the labels strong and fragile were introduced to refer to 
the ability of a liquid in withstanding changes in MRO 
with temperature.E3 In context of bond-ordering model, 
these labels will be referring to the ability of a liquid to 
withstand changes with temperature in bond magnetiza- 
tion mi,. It is worthwhile to mention that in this scheme 
an ideal glass may be characterized as being maximally 
bond ordered which should also imply the least possible 
energy. 



V. SIMULATION: 
PROCEDURE AND RESULTS 

Clock models have been investigated in some detaiL 
both analytically and through numerical methods. tlTO 
Previous Monte Carlo works, have examined the behav- 
ior of 2D clock systems with various values of the pa- 
rameter q. The important results are that for q < 4 the 
system exhibits one second-order transitiop-but for q > 4 
two Kosterlitz-Thouless (KT) transitions^! are present. 
The upper transition temperature is believed to have a 
value approximately equal to the_KT transition temper- 
ature for the continuous modeled T c = 0.89 J/kg- As q 



increases, the lower transition temperature (cxl/q 2 ) ap- 
proaches zero, leaving just one KT transition for the 2D 
XY-model. 

For model Hamiltonian we consider Eq. (6) , which in 
fact is a bond-diluted version of g=8 state clock model, 
involving n.n. interaction with antibonding electronic 
state. For our purpose, the eight possible orientations 
must be interpreted as the possible displacements of an 
atom from its designated equilibrium lattice site, with 
every atom's displacement being the same if the system 
were in a ordered configuration. 

The standard MC importance-sampling method was 
used to simulate the behavior of the system on 
LxL square lattices with periodic boundary condi- 
tions. We performed simulations on lattices of size 
L = 12, 20, 32, and 50. Preliminary work was carried 
out on systems of size L = 12. In every case the temper- 
ature was lowered in steps of 0.05 J/fcs, starting with 
initial value T= 2.00 J/ks- The system was initialized 
in random configuration suitable for high temperatures 
where it is known to be disordered. At every temperature 
system was allowed to equilibrate through a few hundred 
Monte Carlo steps per site (MCS). The data points were 
then acquired by averaging over 40,000 MCS. The data 
were accumulated in several bins, and binned averages 
were used to obtain error estimates for the calculated 
means and also to monitor the state of equilibrium. For 
most cases the estimated statistical errors were less than 
2% of the calculated mean values. In order to test our 
code, we simulated q — 6 state clock model and com- 
pared our results, .ssiith the extensive literature available 
on the subjectO'ta To our satisfaction the agreements 
were impressive. As a note on the calculation of internal 
energy, the (possible) bonds were sorted out in ascend- 
ing order of energy at every time step and were given 
bonding-electrons in that order. 

Many of the quantities calculated have been already 
discussed and given explicit mathematical definitions in 
Sec. IV. In addition, specific heat was obtained from fluc- 
tuations in the total energy: 



C = ((H 2 ) - (H) 2 ) /2Nk B T 2 



(16) 



The bond-related quantities; bond order parameter, 
bond susceptibility, internal energy and specific heat are 
normalized to the number of n.n. bonds 2N, for a system 
of size N = L 2 with periodic boundary conditions. On 
the other hand the structural properties, i.e. structural 
magnetization and susceptibility are normalized to the 
system size N. Results correspond to a hole concentra- 
tion of 40%. 

Perhaps the most interesting feature of the results is 
the conspicuous peak in bond susceptibility, shown in 
Fig. 2, which is an indication of bond ordering, uncorre- 
lated with any long-range structural ordering as is ev- 
ident from structural magnetization and susceptibility 
seen in Figs. 3, and 4. Finite size effects are quite evident. 
The temperature at which bond susceptibility reaches 
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its maximum value is estimated to be 0.9 J/ks, con- 
sistent with the KT transition temperature. At exactly 
the temperature where bond susceptibility maximum oc- 
curs, one observes a maximum in specific heat shown in 
Fig. 5, which should be attributed to large energy fluctu- 
ations associated with the ordering of bonds. Indeed, ex- 
perimental measurements of heat capacity for covalently 
bonded fragile systems such as AS2S3 and B2O3, exhibit 
similar peaks at glass transition which therefore suggest 
bond ordering nature for the glass transition. The in- 
ternal energy for various system sizes is shown in Fig. 6, 
displaying steep slope in bond ordering region. 

Fig. 7 contains the variation of bond magnetization 
with temperature. In bond ordering region, bond mag- 
netization increases rapidly with decreasing temperature 
in spite of the fact that structural magnetization stays 
fairly constant there. This behavior testifies to the earlier 
assertion that a system can undergo substantial bond or- 
dering and hence largely reduce its internal energy with- 
out undergoing any significant structural ordering. 

Unlike the behavior expected from g-state clock model 
in 2D, structural susceptibility (Fig. 4) does not exhibit 
singular behavior at intermediate temperature range 
hence ruling out the possibility of long-range cooperative 
structural ordering. This behavior of structural suscep- 
tibility may be seen in view of the large concentration of 
holes or bond dilution. 



VI. SUMMARY 

We have investigated the existence of order associated 
with bonds, in amorphous systems of interest such as 
vitreous silica. Through our MC simulations, it is found 
that bond ordering may occur independently of struc- 
tural ordering. Bond ordering implies ordering in n.n. 
and n.n.n. distances and thus leads to MRO, which is a 
key aspect of the vitreous state. 

An order parameter for supercooled liquids and glasses 
is introduced. In case of liquid <-> glass transition bond 
magnetization nib may be fruitfully employed in the ther- 
modynamics of glass and any calculation that may in- 
volve coarse-grained parameters. In addition, a new 
identification for glass transition temperature is afforded 
through the variation with temperature of the bond sus- 
ceptibility. 

As a last remark, the sharp rise in the viscosity of glass=. 
forming liquids when cooled toward glass transition,E2l 
may be viewed in terms of strengthening of bonds. 
The fragile glass-formers undergo significant bond or- 
dering through the transition and for that reason their 
viscosity rises dramatically, as opposed to the smooth 
Arrhenius-like behavior of viscosity exhibited by strong 
glass- forming liquids. 
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FIG. 1. The distribution of Si-O-Si bond angle, measured 
by Mozzi and Warren.u 
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FIG. 2. Plot of bond susceptibility vs. temperature for FIG - 5 - Specific heat as a function of temperature for three 

three different system sizes; 20 2 (+), 32 2 (x), and 50 2 (*). different system sizes; 20 (+), 32 (x), and 50 (*). 
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FIG. 3. Plot of structural magnetization vs. temperature. 
Each curve is normalized to the corresponding system size N. 
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FIG. 6. Energy per bond as a function of temperature. The 
total energy E is in units of the coupling J. 



2.0 



£ 1.0 - 



0.0 



+ 


20x20 


X 


32x32 


* 


50x50 


* * 




* g 

+ 5 




+ ? 




* 

* * 
* * 

* 

* M 





1.0 



1.5 



2.0 



FIG. 4. Variation of structural susceptibility with temper- 
ature. Clearly, there is no indication of long-range structural 
ordering in the temperature range shown. 
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FIG. 7. Plot of bond magnetization vs. temperature. There 
appears a sharp variation due to bond ordering, uncorrelated 
with any long-range structural ordering. 
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